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Abstract 

We apply the light-front quantization to the Nambu-Jona-Lasinio model with the 
vector interaction, and compute vector meson's mass and light-cone wavefunction in the 
large limit. Following the same procedure as in the previous analyses for scalar and 
pseudo-scalar mesons, we derive the bound-state equations of a qq system in the vector 
channel. We include the lowest order effects of the vector interaction. The resulting 
transverse and longitudinal components of the bound-state equation look different from 
each other. But eventually after imposing an appropriate cutoff, one finds these two are 
identical, giving the same mass and the same (spin-independent) light-cone wavefunction. 
Mass of the vector meson decreases as one increases the strength of the vector interaction. 
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The light-cone (LC) wavef unction of a hadron is one of the most useful quantities for de- 
scribing the hadron structure in terms of its underlying degrees of freedom [1]. In general, it 
contains information about soft dynamics among quarks, antiquarks and gluons, and one can 
compute various scattering processes involving hadrons in initial/final states, by combining it 
with the hard part of the diagrams. Diffractive vector meson production is one of the typical ex- 
amples of such processes [2]. To compute the amplitude of diffractive electro/photoproduction 
of a vector meson for a wide range of kinematics, one needs to know the LC wavef unction of a 
vector meson with non-perturbative information. In this Letter, we are going to discuss this LC 
wavefunction of the vector meson in a simple model. As another interesting example, E791 ex- 
periment at Fermilab [3] has recently attempted to determine the LC wavefunction (squared) of 
pions through the diffractive pion dissociation process (dijets production) according to Ref. [4]. 
Although it is argued that determination of the pion LC wavefunction from the experimental 
data is actually quite hard [5], it is still true that one cannot compute the amplitude of this 
process without knowing the pion LC wavefunction. Similar experiments are possible in the di- 
jets production from a virtual photon, where, according to the vector meson dominance model, 
the vector meson contribution forms the hadronic part of the photon wavefunction. 

Perturbative calculations provide us with the so called " asymptotic" form of the LC wave- 
functions. For example, the asymptotic form of the pion LC wavefunction is known as well as 
the vector meson's one [6, 7]. However, non-perturbative study is quite few. Lattice simula- 
tion can compute the first few moments of meson's distribution function, but at present they 
are not sufficient to determine the LC wavefunction itself. Therefore, it is very important to 
develop a non-perturbative technique which allows us to directly obtain the LC wavefunction. 
Clearly, the most straightforward and natural framework is the Hamiltonian formalism in the 
light-front (LF) quantization [8, 1]. Before challenging the problem in the real QCD, one should 
be able to learn much from the analyses of simpler models such as the Nambu-Jona-Lasinio 
(NJL) model. Indeed, this model was recently studied by two of the authors within the LF 
quantization [9, 10] and we follow the same procedure to get the LC wavefunctions of vector 
mesons. As is well known, there is a paradoxical situation in the LF quantization. It has been 
asked how one can describe spontaneous symmetry breaking in a formalism having only a trivial 
Fock vacuum. This was answered in Ref. [9] within the NJL model with component fermions 
^N3h = ^{i^ — mg)^ + \Gi[{^'^Y + {^i^^'^Y]. Based on the analogy with the description 
of spontaneous symmetry breaking in a scalar theory on the LF, they found that, still with 
the trivial Fock vacuum, chiral symmetry breaking is described in such a way that one selects 
an appropriate Hamiltonian depending on the phases of the symmetry. In the NJL model, 
different Hamiltonians are originated from different solutions to the constraint equation, which 
exists only in the LF formalism. The "bad" component of the spinor, ip_ (where ■il)± = ^7^7''^^ 
and ■0+ is called "good"^) is not a dynamical variable and is subject to a constraint equation, 
as we will see below. This "fcrmionic constraint" is a nonlinear equation in the NJL model 
and leads to the "gap equation" for the chiral condensate if one adopts an appropriate cutoff. 
Namely, using the parity invariant cutoff Ip^^l < A, one gets 

M-mo GiNA^ ( ^ , 2An 1 

where M = tjiq — is the dynamical mass of the fermion. When the coupling constant 

Gi = GiN / A'k'^ is larger than the critical value (5^"'*"^^'^ = 1/2, the gap equation has a 

^Our notation is the following: = (x° ±x'^)/V2, — [v^^v^ tV\)^v = {v^ ^v\^v\) and d± = d/dx^. 
We use jj,, V for Lorentz indices of four vectors, z, j for transverse coordinates 1, 2, and a,/3 for spinor indices. 



non-zero solution even in the chiral limit tjIq — > 0. This means that the fermionic constraint 
allows for "symmetric" and "broken" solutions corresponding to those of the gap equation. If 
one selects the "broken" solution, and substituting it to the canonical Hamiltonian, one obtains 
the "broken" Hamiltonian. This governs the dynamics in the broken phase and is completely 
different from the Hamiltonian with the "symmetric" solution. In Ref. [9], the fermionic con- 
straint was solved by using the expansion [Indeed, Eq. (1) is the leading order result], and 
they obtained the Hamiltonian in both symmetric and broken phases. They also solved the 
bound-state equations for the scalar and pseudo-scalar mesons^ and obtained their LC wave- 
functions and masses, as well as the PCAC and GOR relations. 



One can of course apply the same procedure for vector states, but we know that the NJL 
model vCnjl does not allow for a bound state in the vector channel [11, 12]. Vector states start 
to bind if one adds the vector interaction so that the attractive force between a quark and an 
antiquark in the vector channel becomes stronger. Therefore, in this Letter, we include the 
vector interaction minimally by adding the following interaction: 

G2 



(: ^7^^^ :f + (: ^7^75* :)'J . (2) 
This interaction, however, makes the fermionic constraint tremendously complicated: 



7+^/'+(^+^_ -I- '^p-^|J+) - i757"^^/'+(^+i75^- + i'-i'^->^+) 
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7i7+^+ : {V'+7iV'- +-0-71^+} : -71757^-0+ : {^/'+7i75'0- + V'-7i75'0+} :] • (3) 

Here we have followed the same operator ordering as in the previous analysis without the vector 
interaction. After rewriting this equation into a bifocal form, one can solve it in the quantum 
level by using the 1/A^ expansion, which is systematically generated by the Holstein-Primakoff 
technique [13]. It turns out that the leading order equation gives the same gap equation as 
Eq. (1). This is natural because we have taken the normal order^ in the interaction (2). 

In the leading order of the expansion, mesonic states are written as constituent states 
with a (dynamical) quark and a (dynamical) antiquark, as was discussed in Ref. [9] for the 
scalar and pseudo-scalar states. Thus, a generic vector state with total momentum P'* and 
hehcity A can be represented as follows: 

|vector;A,P) = dk+ J d^k± (j)x{x, k±) (4) 

X e^(P, A) {r'^i-k, -P + k) + r^(-P + k, -k) Pl^(fe, P - k)\0), 

where A;"*" and k± are the longitudinal and transverse momentum of the quark {x — k'^/P^), 
(f)x{x, k±) is the spin independent part ^ of the LC wavefunction which should be determined 

^Though there was only one flavor in Ref. [9], these physically correspond to the pion and sigma mesons, 
and also the vector meson to be discussed in the present Letter corresponds to the rho meson. Generalization 
to multi flavors should be straightforward. 

'^Normal order is defined with respect to the Fourier modes of the fermionic field. 

^By definition, (f>\{x,k±) should be independent of A, but we retain A because, as we will see below, the 
bound-state equations look different for different A. 



by the dynamics, e^{P,\) is a polarization vector, and B]^p{p^q) is a bosonic operator which 
was introduced to solve the fermionic constraint [9] and corresponds to quark and antiquark 

creation operators ~ h'pdq in the leading order of the expansion. Spin dependent part of 
the LC wavefunction, F'^ — (r+,r~,r*), is determined by the interpolating field of the vector 
meson: 

r+(p,qr)=f+(p,qr) = ^-l, (5) 
r-(p, q) = T-{q,p) = {7lpl7i?i + Mf_,{p^ + q^y + M'} , (6) 

r^(P, q) = -^7l {liqi + m) , r{p, q) = {^{qi + m) 7l . (7) 

The polarization vector is written in the rest frame of the meson P'* = ("^v/ V^, "^v/ V^, 0_l) 
as (see also [14]) 

."(P. A = ±1) = (o. 0, ii, -1) . e^(P, A = 0) = --L. 0,) . (8) 

As we explained before, after solving the fermionic constraint (3) with the nontrivial solution 
M to the gap equation (1) and substituting its solution into the canonical Hamiltonian, 
one gets the "broken" Hamiltonian i^LF- Its exphcit form and the detailed derivation of the 
Hamiltonian are very complicated and will be reported elsewhere [15]. Here only the final 
results are shown. The Hamiltonian is derived order by order of the 1/N expansion, 

OO / 1 \ " 

n=0 



and /i^^-* turns out to be the lowest non-trivial Hamiltonian since is just a constant and 
h^^^ — 0. Therefore, keeping this non-trivial order, one can write the eigenvalue equation for a 
vector state as 

/i(') I vector; P) = ""^^/^ | vector; P) . (9) 

Solving this equation yields both the spin independent part of the LC wavefunction (j)x{x, k±) 
and mass of the vector meson my simultaneously. 

Before going into details, let us briefly discuss the qq states to clarify the procedure we 
perform. In general, the LF energy of the two body state (4) may be schematically written as 

2k+ + 2(P+-k+) +^^^'^^' ^^"^ 

where the flrst two terms are the "kinetic" energies of the quark and the antiquark, and V is 
the potential which allows for a bound state. This form of the energy leads to the following 
bound-state equation: 



' kl + M^ 



x{l — x) 



where we have chosen the vector meson's rest frame, P = (P'^,P^) = (my/v^, 0_l) for sim- 
plicity, and redefined V with some factors included. In the following, since we are interested 



in seeing how the vector interaction (2) affects the vector channel, we will derive the potential 
V up to the leading order of the vector interaction 6*2. We will see that in this leading order 
the potential term V{x, k±; y,p±) is separable with respect to the internal {y,p±) and external 
{x, k±) variables, and actually depends only on y and p±. 

Now let us explicitly show the bound-state equations of the transverse and longitudinal 
components derived from the leading nontrivial Hamiltonian h^'^\ First, for a transversely 
polarized vector meson, a lengthy calculation yields the following potential Vt i^i^) — x/\x\) 
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(12) 



Notice that this is already independent of the external variables x, kj_. Taking the leading con- 
tribution of G2, one arrives at an equation for the LC wavefunction ^^(a;, k±) = 4>x=±i{x, k±): 



rrir. 



kl + M^' 
x{l — x) 
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2 pI + M^- 2y{l - y)p'\ 
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0t(i/>P±)- (13) 



Next, the longitudinal component is much more involved. A longer, but straightforward calcu- 
lation leads to a more complicated potential Vi,: 



X 



G2N 
(27r)3 

2 + 



1-2 



GoN /-i 



(27r)3 
4 {kl + M2 



dz j d'^q_L 



rriT 



kl 



mlx{l -x)- {kl_ + M2) 



G2N 
(27r)3 



x{l — x) 
j dz j d^Qi 



ml + 



nil 



kl 



M 



2 



x{l 



X] 



pi 



M2 



.(14) 



Again, taking the leading term with respect to G2 after careful modification of the bound- 
state equation, we eventually obtain the following simpler equation for the longitudinal mode 

4>h{x, k±) = (px=oix, A';j_): 



A;2 +M2 



^^A^ ^2 4(pi + M2) 
I = - J dyd pi_ y^-^_y^ My,P±) ■ (15) 



x(l 



X) 



It is evident that the right-hand side is again independent of the variables x, k±. Since the step 
from Eq. (14) to Eq. (15) is a bit nontrivial, let us show the easiest way to derive Eq. (15), 
which is however less systematic. First of all, if one ignores the G2 dependent term in the 
second line of Eq. (14) that gives the higher order in G2 and thus can be ignored anyway, one 
immediately finds 
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Then one integrates the resulting bound-state equation over x and A;^, obtaining the following: 



/ 



dxd k\ < mr 
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One can modify the bound-state equation by using this integral in the right-hand side of it. 
Finally, taking the leading term with respect to G2, one obtains Eq. (15). It should be noted 



that the above equation is consistent with Eq. (15) since it is the integration of Eq. (15) over 
X and k^. 

At first glance, the above two eigenvalue equations (13) and (15) look different and thus 
seem to give different masses for the transverse and longitudinal vector mesons. This is of 
course physically unacceptable, and as we will verify soon, these equations are essentially the 
same and give the same mass — m-^. This equivalence will be achieved after one specifies 
cutoff scheme. It is not hard to identify the origin of this (fake) difference with the lack of 
Lorentz covariance in the LF formalism. 

Nevertheless, even at this stage, one can see that solutions to the above equations are 
an identical function of x and k^. To this end, it should be noted again that these two 
equations have a very simple structure: Their right-hand sides depend on neither x nor k^. 
This immediately implies that the solutions should be 

"^T/L x{\-x) 

Thus, if the masses of transverse and longitudinal vector mesons coincide with each other, then 

so do the LC wavefunctions (Ct/l are determined by the normalization). The LC wavefunction 
(16) has a peak at x = 1/2 as we will see explicitly below. Since our description of the vector 
states is with respect to the quark (antiquark) having a dynamical mass M 7^ 0, this shape of 
the LC wavefunction implies the constituent picture. 

Now let us verify that both the two equations (13) and (15) derive the same equation for a 
vector meson mass my. Inserting the solution (16) into these equations, one arrives at integral 
equations for my. We evaluate the integrals by introducing the "extended parity invariant 
cutoff' [9] which is actually equivalent to the Lepage-Brodsky cutoff [6]: 

Indeed, this is a natural extension of the parity invariant cutoff in the two body sector, K'^K~ < 
where are the sum of (on-shell) quark and antiquark longitudinal momenta and energies 
[K+ =p+ + (P+ - p+) = P+, K- = (pI + M^)/2p+ + (pI + M2)/2(P+ - p+)]. This cutoff 
apparently preserves transverse rotation and parity symmetry separately, but in fact it does 
work better. First, it also respects the usual three dimensional space rotation [16]. Thus 
one can relate the above cutoff A to the 3-momentum cutoff J2i=i,2,3{P^T < -^Im through 
2(A|]v[ + M^) = A^. Next, the cutoff K^K~ < is invariant under the boost transformation 
K"^ e^^K^, which is necessary for the relativistic formulation. 

In Ref. [9], the parity invariant cutoff was specified as \K^\ < A which contains two inde- 
pendent conditions. In actual calculations, however, the authors of Ref. [9] utilized only the 
Lepage-Brodsky cutoff which is obtained by conbining the two conditions^. Namely, the con- 
dition IK^^I < A was introduced only to derive the Lepage-Brodsky cutoff. However, putting 
the cutoff on the longitudinal momentum < A for the total momentum is not preferable 
from the viewpoint of boost symmetry. Thus, in the present paper, we redefined the parity 
invariant cutoff in the two body sector by K^K^ < A^. On the other hand, there is no prob- 
lem in putting |j9''=| < A in the gap equation (1) because the momentum is not the external 

^Therefore, body of the calculations in Ref. [9] is correct, while the derivation of the Lepage-Brodsky cutoff 
was not appropriate. 



momentum but the internal one to be integrated out. Indeed, the gap equation comes from the 
zero longitudinal momentum of the fermionic constraint written in the bilocal form. 
Now the integral in the equations is replaced as follows^°: 

jdyj d^p^ dy ndipl) , (18) 

with 



1±3 2M2 ^ ^ 

Then one can exphcitly prove that the two equations from Eqs. (13) and (15) do give the same 
equation that determines the mass ruy = nir^ = m^. It is very important to recognize that 
we can derive this equation simply by inserting the LC wavefunction (16) into the bound-state 
equations with the above cutoff. We have just evaluated the right-hand sides of the equations. 
Since the LC wavefunction (16) is a direct consequence of the bound-state equations, to obtain 
the same equation for my means that the original equations are also equivalent to each other. 
The explicit form of the equation for the vector meson mass is given by 



1 _ 2 
G^" 3 




arctan ■ 



(20) 



where we have defined a dimensionless couphng constant G2 — G2N / At:'^ and r is (square 
of) the ratio of the vector meson mass to the threshold mass 2M: 



my ^ ^ 

2M 



(21) 



A physical bound-state appears only when the ratio r is in the range < r < 1. Equation 
(20) has a solution in this region when the strength of the coupling constant G2 is in the range 
Gt"""^ < G2 < &r'''^ defined by 

Gr^-^{/5 + (l-/?^)ln(i^)p, Cr-)^^-^. (22) 

Two hmiting cases G2 — ^2™°^ and G2 — G2^^^ correspond to r = 1 (loose binding limit), and 
r = (tight binding limit), respectively. When M/A — > (/3 — > 1), the physical bound-state 
region shrinks G2^^\ while it becomes wider as M/A grows large. The existence of 

^(min) consistent with the observation that there is no bound state in the NJL model without 
the vector interaction. Similar behaviors have been found in Ref. [17]. 

In Figure 1, a numerical solution to Eq. (20) is shown as a function of G2, where the 
constituent quark mass is taken to be M/A = 0.4 (/3 = 0.82) as an example. As we expect, the 
bound state appears for G2 larger than some critical value and the mass starts to decrease from 
the threshold value 2M as one increases the strength of the vector interaction. The value of 
critical coupling constants are exactly the same as the values predicted by analytic calculation. 
When (3 = 0.82, they are G^'''^ = 0.95, G^""^"^ = 2.72. 



^'^Alternatively, one can regard that the LC wavefunction (16) has support defined by the Lepage-Brodsky 
cutoff. 
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Figure 1: my/2M as a function of G2- Constituent quark mass is M/A = 0.4. A bound state 
appears only in the regime Ga™™^ = 0.95 < G2 < G^^^"^^ = 2.72. 



It is also interesting to see my as a function of Gi and 6*2 • The dependence on Gi enters 
only through the constituent mass M [see Eq. (1)]. In Figure 2, we show the vector meson mass 



m 



V 



as a function of Gi = GiNA"^ /Att'^ and G2 in the broken phase Gi > Gi 



(critical) 



1/2 and 



in the chiral limit. As the coupling constant Gi becomes large, the constituent mass (namely, 
the chiral condensate) becomes large. For fixed cutoff A, this means to increase the value M/A 
and thus enlarges the bound-state region. This can be seen clearly in the figure. 

In Figure 3, we compare the LC wavefunction of the vector meson (16) with that of the 
pseudo-scalar meson (pps{x,k±) oc (m|g — {k^ + M^)/x(l — x))^^ which is the result of the 
previous analysis where the vector interaction was not included [9]. The transverse momen- 
tum k± is set to be zero for simplicity and the wavefunctions are normalized at x = 1/2 for 
comparison^^. We have chosen the masses to be (mps/2M)^ = 0.01 and (mv/2M)^ = 0.95 as 
a typical case in the chirally broken phase with non-zero current quark mass mg 7^ 0. In spite 
of the absence of the vector interaction in the previous analysis, this comparison makes sense 
because the effects of the vector interaction on the mass of a pseudo-scalar meson is small. 
Indeed, after the Fierz transformation, the vector interaction will generate terms proportional 
to the original Gi interaction, but of the higher order in 1/A^. Notice that the LC wavefunc- 
tion of the pseudo-scalar meson has the same functional form of x and k± as that of Eq. (16). 
Therefore, difference of the shape is due to different values of the mass. The constituent picture 
works better in the vector meson than in the pseudo-scalar meson. 

Before concluding the paper, let us show one more evidence for the equivalence between 
the transverse and longitudinal equations. As we already mentioned, the superficial difference 
comes from the fact that the Lorentz covariance (in particular, the 3 dimensional rotational 
invariance) is not manifest in the LF coordinates. If one works in a framework with obvious 
rotational invariance, then there should be no difference between transverse and longitudinal 
components. Indeed, in the leading order of the 1/A^ expansion, it is possible to derive the 
same equations from the covariant Bethe-Salpeter (BS) equation: 



IBS 



G2N 



d^k tr 



7i0^^(fc;P) 



^- M + ie 
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^-f -M 



(23) 



-'^-'^ If one includes the Lepage-Brodsky cutoff into the definition of the LC wavefunction, support of x for 
fcx = is given by y_ < a; < j/+ with y± given by eq. (19). 




comparison. 



where (j)^^{q; P) is the amputated BS amphtude for the vector channel (P and q are the total 
momentum and relative momentum of a quark and an antiquark system, respectively). The 
LC wavefunction is obtained by the LF energy integral of the BS amplitude. Here, however, 
wc have to be much careful because we must specify the cutoff scheme to make the equations 
well-defined. When we impose the cutoff, it is of course (and again) desirable to maintain the 
symmetries such as the rotational invariance. For example, in Eq. (23), we can use the so called 
3-momentum cutoff scheme which respects 3-dimensional rotational invariance. Then, it turns 
out that the resulting two equations are equivalent to Eqs. (13) and (15) respectively, with the 
parity invariant cutoff (17) if one uses the following variable transformation y <-> {Ek + kz)/2Ek, 

with Ek = \Jk'^ + M"^ and k b cing the three momentum here. Also, one needs to replace the 3- 
momentum cutoff A3M by the extended parity invariant cutoff A using the relation shown before. 

To summarize, we have applied the LF quantization to the NJL model with the vector 
interaction, and obtained the eigenvalue equations for vector meson's LC wavcfunctions. Due 
to the addition of the vector interaction, the vector state becomes a bound state. At first glance, 
transverse and longitudinal components of the bound-state equations look different from each 
other, but eventually after imposing an appropriate cutoff scheme, one finds these two coincide 
with each other. Mass of the vector meson decreases as one increases the strength of the vector 
interaction. This behavior is consistent with Rcfs.[17, 18] which also treated the vector meson 
within the NJL model with the vector interaction. 

Once we obtain the LC wavcfunctions, we can compute various physical quantities. One 
of such important quantities is the physical form factors. For the pseudo-scalar case, this was 
done [19] in the NJL model with two flavor quarks. Similar analysis can be done in the vector 
meson case and will be reported in the future publications [15]. 

Our non-perturbative approach plays a complementary role to the perturbative calculation 
of the asymptotic form [7], because our LC wavefunction is expected to describe that of low en- 
ergy scale (Of course we have to include flavor degrees of freedom, which is straightforward. See 
for example, [16]). It is interesting to flnd a way to interpolate these two different approaches. 
One of the possible ways for this problem is to include the effects of gluon propagation between 
a quark and an antiquark. This is easily incorporated by using the nonlocal current current in- 
teraction j'^Dni.j" instead of the point interaction in the NJL model. Such kind of interpolation 
between low and high momentum regimes is used in various situations. It is also interesting to 
apply our approach to the heavy quark system such as J/V' or T. The use of the nonlocal inter- 
action is also convenient from the technical point of view. In the presence of the gauge held, one 
needs to consider the longitudinal zero modes of the gauge fleld. In spite of the favorable aspect 
that the nontrivial vacuum structure such as the theta vacua may be attributed to the dynamics 
of gauge zero modes, inclusion of the gauge zero modes makes the canonical structure terri- 
bly complicated [20]. Thus, for the problem of quarkoniums where one can ignore the vacuum 
physics, it is easier to replace the gauge fleld effects by the non-local current-current interaction. 
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